CONJUGATE POINTS IN NILPOTENT SUB-RIEMANNIAN PROBLEM 
ON THE ENGEL GROUP 
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Abstract. The left-invariant sub-Riemannian problem on the Engel group is considered. This prob- 
lem is very important as nilpotent approximation of nonholonomic systems in four-dimensional space 
with two-dimensional control, for instance of a system which describes movement of mobile trailer 
robot. We study the local optimality of extremal trajectories and estimate conjugate time in this 
article. 



1. Introduction 

This work deals with the nilpotent sub-Riemannian problem on the Engel group with growth vector 
(2,3,4). Four-dimensional optimal control problem with two-dimensional control is stated as follows: 
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Since the problem is invariant under left shifts on the Engel group, we can assume that the initial 
point is identity of the group qo = (xo,yo, zq,vo) = (0, 0, 0, 0). 

The paper continues the study of this problem started in the work [5]. The main result of [5] is 
upper bound of the cut time (i. e., the time of loss of global optimality) along extremal trajectories 
of the problem. The aim of this paper is to investigate the first conjugate time (i. e., the time of loss 
of local optimality) along the trajectories. We show that the function that gives the upper bound of 
the cut time provides the lower bound of the first conjugate time. In order to state this main result 
exactly, we recall necessary facts from the previous work [5]. 

Existence of optimal solutions of problem {[TJ) ([3j) is implied by Filippov's theorem [3]. By Cauchy- 
Schwarz inequality, it follows that sub-Riemannian length minimization problem ([3]) is equivalent to 
action minimization problem: 
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Pontryagin's maximum principle [HJ [6] was applied to the resulting optimal control problem dTJ) , (EJ) , 
Abnormal extremals were parameterized. 
Denote vector fields at the controls in the right-hand side of system ([!]): 



(i,o,-|,of , 
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and the corresponding linear on fibers of the cotangent bundle T*M Hamiltonians /ij(A) = (X,Xi(q)), 
A G T*M, i = 1,2. Normal extremals satisfy the Hamiltonian system 

A = H(X), A G T*M, (5) 

where H = \ (h\ + ftf) • 

The normal Hamiltonian system ([5]) is given, in certain natural coordinates, as follows on a level 
surface {A G T*M \ H = ±}: 

9 = c, 6 = —as'm9, a = 0, (6) 
g = cos6>Xl(<?) + sin^X 2 (g), g(0) = g - 

The family of all normal extremals is parameterized by points of the phase cylinder of pendulum 
C = |a G T* q M I H(X) = i| = {(#,c,a) | G S\ c,a G M} , 

and is given by the exponential mapping 

Exp : N = C x M.+ — >■ M, 
Exp(A,t) = q t = (x t ,y t ,z t ,v t ). 

Energy integral of pendulum ([6]) is expressed by E = ^ — a cos 9. The cylinder C has the following 
stratification corresponding to the particular type of trajectories of the pendulum: 

C = U 7 i=1 d, a n Cj = 0, i + j, A = (0, c, a), 
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Extremal trajectories were parameterized by elliptic Jacobi's functions for any A G C in the paper [5]. 
This parameterization was obtained in natural coordinates ((p,k,a), which rectify the equations of 
pendulum: tp = 1, k = 0, a = 0. 

Further, in the work [5] discrete symmetries of the exponential mapping were described. The 
corresponding Maxwell sets were constructed. On this basis was obtained the main result of the 
paper [5J, Theorem HJ which gives upper bound of the cut time along extremal curves 

^cut(A) = sup{i > | Exp(A,s) is optimal for s G [0,i]}. 

Define the following function £m AX : C — > (0, +oo]: 

A G Ci => 4ax = min(2pi, 4K)/<r, 
A G C 2 => = 2Kk/a, 

A G C 6 => tJ^Ax = — , 

A G C 3 U C 4 U C 5 U C 7 =>- 4iAX = +oo. 
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/dt 
— ; p\{k) G (K(k), 3K(k)) is the first positive root of the 
VI - k 2 sin 2 t 
o v 

function f z (p,k) = dnp snp + (p — 2E(p))cnp; dnp, snp, cnp are Jacobi's functions [14J; E(p) = 
/ P dn 2 tdt. 

Theorem 1 ( [5], Theorem 3). For any A G C 

tcut(A) < 4ax(A). (7) 

We study the local optimality of extremal trajectories and estimate conjugate time in this article. A 
point q t = Exp(A, t) is called a conjugate point for qo if v = (A, t) is a critical point of the exponential 
mapping and that is why qt is the corresponding critical value: 

d u Exp : T U N — > T qt M is degenerate, 

i. e., 

d(x,y,z,v) , = 
d(6,c,a,t) {U> 

Note that t is called a conjugate time along extremal trajectory q s = Exp(A, s), s > 0. 

Here and below we denote by — — - — - — ^— the Jacobian of the exponential map 

o(p, c, a, t) 
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Due to the strong Legendre condition, for any normal extremal there exists a countable family of 
conjugate points. Besides, conjugate times are separated from each other (see Section [3]). The first 
conjugate time along the trajectory Exp(A, s) is denoted by 

*conj = m i n {* > | t is a conjugate time along Exp(A, s), s > 0} . 

The trajectory Exp(A,s) loses local optimality at the moment t = £<!: on j(A) (see [3J). Our main aim 
is to prove the following lower bound of the first conjugate time. 

Theorem 2. For any A G C 

*conj(^) > *MAx(^)- ( 8 ) 

In Sections HH3 we prove the inequality (jHJ) , A G Cj for alH = 1 , . . . , 7 (see Theorems [H EJ O [7]) . 

2. Conjugate time and symmetries of the exponential mapping 

Normal Hamiltonian system for the considered problem has the following symmetries (see [5]): 
reflection 

(9, c, a, x, y, z, v, t) i->- (6 — n, c, —a, —x, —y, z, —v, t) (9) 

and dilations 

g Q/ 3 

(9,c,a,x,y,z,v,t) ^ (e,—,-,^x,y^y,-fZ,^v,y^yt), 7 > 0. (10) 

V7 7 

We consider the corresponding symmetries of the exponential mapping and their action on conjugate 
points. 
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2.1. Reflection. Define the action of reflection in preimage and image of the exponential mapping 
according to ([9]): 

i : N — > N, i{v) = i(9, c,a,t) = v = (0 — it, c, —a, t), 
i : M M, i(q) = i(x,y,z,v) =q = (-x, -y, z, -v). 

Existence of symmetry ([9]) of Hamiltonian system implies that the reflection i is the symmetry of the 
exponential mapping: Expoi = % o Exp. (It is easily shown that i and the reflection e 4 coincide [5]). 
Hence we obtain dExpodi = diodExp. The reflection i is non-degenerate (Ker di = 0) and therefore 
v = (A, t) = (9, c, a, t) is a critical point of Exp if and only if v = i{y) = (A, t) = {6 — tt, c, —a, t) is a 
critical point of Exp. And so iJ onj -(A) = tl on ^(X). Using the definition of Maxwell time t^x ( see 0> P' 
7.6.), we get easily similar equality %Ax(^) = ^max(^)- Therefore it is enough to prove the necessary 
inequality (JH]) can be proved only for a > 0. 

2.2. Dilations. According to formula (|10p define the action of dilations in preimage and image of 
the exponential mapping: 

~ ^ / c a \ 

$ 7 :M^M, $ 7 (g) = * 7 (x, y, 2, «) = q = (y/yx, yfiy, 7*j 7M i 7 > 0. 

These formulas define the action of multiplicative Lie group M+ in iV and M, s. t. 

Exp o<]? 7 = $ 7 o Exp V7 > 0. 

Thus, there is a one-dimensional symmetry group G = {3> 7 |7 > 0} of the exponential map. 

It is easy to see that the symmetries preserve the sets of critical points and critical values of the 
exponential mapping. 

Lemma 1. 1) If q E M is the critical value of Exp corresponding to a critical point v £ N then 
&y(q) is also the critical value o/Exp corresponding to the critical point < & 7 (^) for any 7 > 0. 

2) Lety>0, \ = (0,c,a),\=(9,^,^ EC. Then t^(A) = ^4njP0 ■ 

Proof. 1) follows from the equality dExpo(i<l> 7 = d<]? 7 o a! Exp. 

2) follows from 1). □ 

Let a > 0. Suppose 7 = a; then from Lemma CO we get the following: 

1 1 1 ( c 

^conj (^J ^1 ^0 /^^conj ( @ i r^i 1 



From the definition of Maxwell time ijyrAX' a similar equation follows: 

4ax(^i c i «) = "^*MAX (#i "7=) 1 ) • 

V« V v a / 

Therefore, it is sufficient to prove the required inequality ([8]) in two cases: for a = 1 and a = 0. 

2.3. Transformation of Jacobian of the exponential mapping. For a fixed A = (9,c,a), 

<9fx, y^ vj 

conjugate times are roots t > of the Jacobian ' — - — First, we transform this Jacobian by 

o(v,c,a,t) 

using the symmetry group G = {3> 7 |7 > 0}. The coordinate expression of the equation Exp o$ 7 (A, t) = 
<1> 7 o Exp(A, t) is 
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Differentiating this equation w.r.t. 7 for 7 = a = 1, we get 

cdq dq t dq 
~2~dc ~ da + 2~di 

Therefore, when a = 1 

d(x,y,z,v) 



x y 3 
— ,—,z,-v 
2' 2' '2 
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d(9, c, a, t) 



j , dq dq dq dq 
dO' dc' da'' dt 



dq dq dq 
det, M'9?9? L 
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Conjugate points and homotopy 

In this section we recall some necessary facts from the theory of conjugate points in optimal control 
problems. For details see, e.g., [HO [13]. 

Consider an optimal control problem of the form 

q = f(q, u), q£M, u 6 U C R m , (11) 

9(0) = qo, q(ti) = qi, ti fixed, (12) 

h 

J = J <p(q(t),u(t))dt min, (13) 



where M is a finite-dimensional analytic manifold, f(q, u) and <p(q, u) are respectively analytic in (q, u) 
families of vector fields and functions on M depending on the control parameter u G U, and U an 
open subset of M m . Admissible controls are u(-) € LqoQO, ti], U), and admissible trajectories q(-) are 
Lipschitzian. Let 

huW = (A, f(q, u)) - <p(q, u), A G T*M, q = 7r(A) G M, 

be the normal Hamiltonian of PMP for problem (]lip -( fT3l) . Fix a triple (u(t), Xt, q(t)) consisting of a 
normal extremal control u(t), the corresponding extremal Aj, and the extremal trajectory q(t) for the 
problem (pT]) - f[T3]) . 

Let the following hypotheses hold: 

d 2 h u 

(HI) For all X £ T*M and u G U , the quadratic form — — ^(A) is negative definite. 
(H2)For any X G T*M , the function u 1— > h u (X), u G U , has a maximum point u{X) G U : 



max h u (X), 



X G T*M. 



(H3) The extremal control u(-) is a corank one critical point of the endpoint mapping. 

(H4)All trajectories of the Hamiltonian vector field H(X), X G T*M, are continued for t G [0, +00). 

An instant > is called a conjugate time (for the initial instant t = 0) along the extremal A( if 
the restriction of the second variation of the endpoint mapping to the kernel of its first variation is 
degenerate, see [3] for details. In this case the point g(i*) = tt(A^) is called conjugate for the initial 
point go along the extremal trajectory q(-). 

Under hypotheses (H1)-(H4), we have the following: 
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(1) Normal extremal trajectories lose their local optimality (both strong and weak) at the first 
conjugate point, see [3]. 

(2) An instant t > is a conjugate time iff the exponential mapping Exp t = 7r o e tH is degenerate, 
see pQ. 

(3) Along each normal extremal trajectory, conjugate times are isolated one from another, see |13j . 
We will apply the following statement for the proof of absence of conjugate points via homotopy. 

Theorem 3 (Corollary 2.2 [7]). Let (u s (t),\f), t G [0, +oo), s G [0, 1], be a continuous in parameter s 
family of normal extremal pairs in the optimal control problem (jlip -(|13 p satisfying hypotheses (HI)- 
(H4). 

Let s i— > if be a continuous function, s G [0,1], if G (0,+oo). Assume that for any s G [0,1] the 
instant t = if is not a conjugate time along the extremal Af . 

// the extremal trajectory q°(t) = vr(A^), i G (0, t®], does not contain conjugate points, then the 
extremal trajectory q ,1 (i) = vr(A|), i G (0, t\], also does not contain conjugate points. 

One easily checks that the sub-Riemannian problem (JTJ) , ([2J) , satisfies all hypotheses (H1)-(H4), 
so the results cited in this section are applicable to this problem. 

4. Estimate of conjugate time for A G C\ 

4.1. Evaluation of Jacobian. We use the elliptic coordinates (ip, k, a) in C±, see [5]. For a fixed 
A = (9, c, a) G Ci, conjugate times are roots i > of the Jacobian 

d(x,y,z,v) 



J 



d(t,(p,k,a) 



d(x, y, z, v 

We transform this Jacobian in the same way as the determinant n/ J — ■ — — 4- in Subsection 



get 



d(9, c, a, t) 



then 



J 



dx dx dx 

dt dip dk 

dy dy dy 

dt dip dk 

dz dz dz 

dt dip dk 

dv dv dv 

dt dip dk 
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(Here and below we assume a = 1 according to Subsec. 12.21 ) Explicit calculation of the function 
with the use of the parameterization of the exponential mapping obtained in [5] gives the following 
expression of the determinant: 



J 
R 



R- Ji 



32 



k(l — k 2 )(l — k 2 sin 2 u\ sin 2 U2) 2 
J\ = do + c?2 sin 2 U2 + d& sin ui , 
d i = d i (u 1 ,k), i = 0,2,4, 
u\ = am(p, k), U2 = am(r, k), 



7^0, 



P = 

d 
di 



t 

2' 

= a\ ■ sinui 
= k 2 a2 ■ f zu 



t 



do + d 2 = -a 2 ■ f z 



(14) 

(15) 

(16) 
(17) 
(18) 
(19) 
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ai = - |4(1 - k 2 ) cosni (1 - 2k 2 sin 2 ui) \J\-k 2 sin 2 niF 2 (ui) + +4/c 2 costti sin 2 tii (1- 

3 

-k 2 sin 2 tii) 5 +4sinni (l - /c 2 sin 2 u x ) £ 3 (tti) - 2 (l - /c 2 ) sinui (l - k 2 sin 2 m) F 3 (tii) + 
+ 2F(«i) (sinm - 2/c 2 (3 - 2k 2 ) sin 3 Ul + k A (5 - Ak 2 ) sin 5 Ul ) ++F 2 (-u 1 )(2 (4/c 2 - 5) (1- 

— A; 2 sin 2 uij sinuiF(ui) + 6 cos tii (1 - 2k 2 sin 2 m) \A - k 2 sin 2 tii) + E( Ul ) (2 (4/c 2 - 
—5) cos u\ (l — 2A: 2 sin 2 u\ ) a/1 - A; 2 sin 2 u x F(u x ) + 8 (l - k 2 ) (l - k 2 sin 2 tii) sinniF 2 (tii) — 

-2 (l + A; 2 + 3/c 2 cos(2«i)) sintti (l - k 2 sin 2 m) 

02 = - cos ui ((F(ni) - F(ni)) 2 + k 2 F( Ul ) (2E(u x ) - F(u x ))\ - 

- sintii\/l - k 2 sin 2 tii - (l - & 2 ) ^(«l)) , 

/zu = sin u\\/ 1 — k 2 sin 2 tti + (F(u\) — 2E(u\)) costti. 

Denote 

x = sin 2 tt2- (20) 

4.2. Conjugate points as k — > 0. We show that extremals corresponding to sufficiently small 
values of the parameter k have no conjugate points for t < ^max(^)- 
The function J\ has the following asymptotics as k — >• 0: 

J\(u\, x, k) = k 2 Ji (tti, s) + o(k 2 ), x = sin 2 tt2, 
J° (tii,x) = d%(ui) + a%(u x )x, 

c^(tii) = - sintii(2tii sin tii + 3u 2 costti + ui sin 3 tti — 6tti sintii + 3costii sin 2 tti), 

^2(^1) = ~ u i sin 4 ui — 2 cos tii sin 3 tti + 3ui sin 2 tti — uf . 

4-2.1. Auxiliary lemmas. We use the following statement to obtain the necessary bounds for func- 
tions. 

Lemma 2. Let real analytic functions f(u),g(u) satisfy on (0,uq) the conditions: 

/(tt)^O, g(u)>0, >0, (21) 

lim M = 0. (22) 
u->0 g[u) 

Then f(u) > for u € (0, u ). 

If functions / and g satisfy conditions (|2ip . (|22p . then we say that 5 is a comparison function for / 
on the interval (0, tto). 

//Y 

Proof. The function — is real analytic, thus it either has isolated zeros or is identically zero. It 
\9j 

( f\ f f 

is not hard to prove that the second case is impossible: if I — =0 then — = const, hence — = 

\9 J 9 9 

f '(^) 

(because lim = 0), whence / = 0, this contradiction proves the case. 
m->o g(u) 
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( f\' f 

So the function I — J > has isolated zeros therefore — strictly increases for u G (0, uq). The 

f C^) 

inequality > follows from the equality ([22]) . so the inequality /(u) > for u G (0,uo) follows 

from g(u) > 0. □ 

We use the following statement to estimate a quadratic polynomial. 

Lemma 3. // /(0, y) > 0, /(l, y) > and a(y) < /or i/ie function f(x, y) = a(y)x 2 + o(y)x + c(y) 
with y G (0,y ), then f(x,y) > for y G (0,y ),z G [0, 1]. 

Proof. We obviously have the inequality f{x,y) > for x = and x = 1, y G (0, yo)- Since a(y) < 0, 
it follows that f(x,y) is convex w. r. t. the variable x (possibly not strictly). Consequently we get 
f(x,y)>0 for ye(0,y ),xe [0,1]. □ 

In the following three lemmas we analyze the sign of the function (ui, x), which is dominant term 
of asymptotics for the function Ji(u\,x,k) as k — > 0. 

Lemma 4. TTie function c?q(ui) = | sin it i (2iii siniti(uf — 3) + 3 cos u\(u\ + sin 2 u\) +u\ sin 3 ui) < 
/or u\ G (0,7r). 

Proof. We show that the function <?(iii) = sinui(sintii — u\ cos tii) is a comparison function for — o]q{u{) 
for u\ G (0, 7r). 

The inequality o![](ni) ^ follows from the expansion d[](ui) = —^f^u 11 + o(u 11 ). 
If u\ G (0,7r), then sintii > 0. Further, 4>{u\) = sinui — Uicosni > for u\ G (0,7r), since 
<^>(0) = 0, 0'(ni) = uisinui > 0. Therefore g{u\) > for u\ G (0,7r). 
Finally we get the equalities 

■d§(«i)V _ (2(-l + u 2 + cos(2iii)) +msin(2m)) 2 



<?(ui) / 4(siniii — ui cos ui) 2 

So g(u\) is a comparison function for — c^(ui) thus it follows from Lemma [2] that c£[j(ui) < for 
niG(0,vr). □ 

Lemma 5. If u\ G (0, ir), then 

do(ui) + d\{u\) = — ( — 2u\ — sinni(— 2u\ sinui + u\ sin 3 u\ + costii(— 3« 2 + sin 2 -ui))) < 0. 

Proof. We check that the function g{u\) = smui(siniti — u\cosui) is a comparison function for 
-(d° ( Ul ) + d° 2 ( Ul )) for Ul G (0,tt). 

The inequality d^iui) + d^{u\) ^ follows from the expansion d^{u\) + d^iux) = —^u 9 + o(u 9 ). 

Note that g(ui) > for m G (0,tt) (see the proof of LemmaS]). 

Also, there hold the equalities 

(dg(ui) + c^(«i))V _ (-2ni + sin(2u 1 )) 2 



g(ni) / 4sin 2 ui 

and 

-(dg(m) + dg(m)) 4 5 5 

^ " 45^ + 

Finally, ^(tii) is a comparison function for — + d^{ui)) therefore it follows from Lemma [2] 

that dg(tti) + ^(iti) < for m G (0, tt). □ 



8 



Lemma 6. For any u\ G (0, ir), x G [0, 1], we have J®(ui,x) < 

Proof. If tii ^ (0, 7r), then Jj(«i, 0) < (see Lemma SJ), J° 1) < (see Lemma [5]), and so it follows 
from Lemma [3] that J° (iti,x) < for u% G (0, ir),x G [0, 1]. □ 

4-2.2. Estimate of conjugate time as k — > 0. 

Proposition 1. There exists k G (0,1), s. t. for any k G (0, k), u% G (0,tt), x G [0,1] we have 
Ji(u%,x, k) < 0. 

Proof. Assume the converse. Then there exist sequences {k n }, {u™}, {x n }, n G N, s. t. k n G (0,1), 
k n — > 0, u r { G (0,7r), x n G [0, 1], and J\ (u^,x n , k n ) > for all n G N. By passing to subsequences, we 
can assume that — > u\ G [0, 7r], x n — )• x G [0, 1]. 

1) Let u\ G (0,7r). From Lemma EJ we get Ji(«i,x) < for all U\ G (0,7r), x G [0,1]. Thus 
Jl(ui, x n , k n ) = k 2 (j° (n™, x n ) + o(l)) < for large values of n, a contradiction. 

2) Let ui = 0. As k 2 + u 2 — > we have 

d ° = "4^5 A;2 ^ 1 + 0(fc2 ^ 1) ' 

d 2 =-4F fc2 «l+°(^l)> 
135 

d4 = ^gfc 4 u? + o(fc 4 «?)- 

2.1) If x / 0, then J\ = — ^k 2 u\x + o{k 2 u\). Therefore Ji(ti™, x ra , k n ) < for large values of n, a 
contradiction. 

2.2) If x = 0, then 

Ji = -J^^l 1 + ^M 1 ) - t^ 2 ^ + °( fc2 ^) 

and Ji(n", x n , fc n ) < as n — )• oo, a contradiction. 

3) Let ui = 7T. As /c 2 + (7r — ni) 2 — ^ we get 

3 

d = --7r 2 /c 2 (7r - ui) + o {k 2 (ir - u{j) , 

d 2 = -TT 3 k 2 + o{k 2 ), 
d 4 = 7r 3 k 4 + o(k 4 ). 

3.1) If x 0, then J\ = k 2 (— 7r 3 x + o(l)) . Whence J\ (ii™, x n , k n ) < as n — > oo, a contradiction. 

3.2) If x = 0, then 

3 

Jl = --7r 2 A; 2 (7r -m) + o (k 2 (it - u{j) - ir 3 k 2 x + o(k 2 x). 

Hence J\(ui,x n ,k n ) < as n — > oo. The contradiction completes the proof. □ 

Going back from the variables (u\,x, k) to (t, ip, k) by formulas (|20"j) . (JHJ), (fl~5j) . we get the following 
statement from Proposition [TJ 

Corollary 1. There exists k G (0,1), s. t. for any k G (0, k), (p G R, an arc o/ i/ie extremal curve 
Exp(A,i), A = (ip,k,a), t G (0, t^ AX (A)) ; does not contain conjugate points. 
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4.3. Conjugate points at t = ^max - ^ n this subsection we find conditions, under which Maxwell 
time t = is a conjugate time. Let us recall that ijviAx(^) = mm (2^(A;), AK{k)) for A G C%, a = 1, 
where p = p\{k) G (if, 3K) is the first positive root of the function /^(p, k) = dnpsnp+(p— 2 E(p)) cnp 
(see [5]). 

It is shown in [llj that 

ke(0,k o ) p\{k) G (3K, 2K), 

k = k => j£(fc) = 2K, 

fe e (/co, l) P l(k)e(K,2K), 

where ko ~ 0.9 is the unique root of the equation 2E(k) — K{k) = 0. Therefore 

UK{k) for fc G (0,fco ] , 
\2j£(fc) for k G [fco, 1) . 

Changing the variable t by u\ = am (|, fc), we get 



*maxM 



u max(^) 



7r for k G (0, feo] j 
u\{k) for fe G [ko, 1) 



where ui = u\{k) = am(pj(fc), fc) G (§,^f-) is the first positive root of the function f zu {u\,k) = 
/ z (am iii, fe)- 

Lemma 7. T/ie function f 3 (k) = E(k) + 2(fc 2 - l)^(/c)if(A;) - (k 2 - l)K(k) > on the interval 
k G (0,1). 

Proof. Let us prove that the function g(k) = 1 — k 2 is a comparison function for /^(/c) on the interval 
k G (0,1). 

The inequality fs(k) ^ follows from the expansion fs(k) = \k 2 + o{k 2 ). 
Notice that g(k) > for k G (0, 1). 

/ 3 (fc)V 2^ 2 (A;) J / 3 (fc) ^ 2 , 2 , n2 , 



Also, we have the equalities — 7— = 7715 -77 and — 7— = — k +0 

V 9{k) J (k 2 - l) 2 g(k) 4 

Finally, g(k) is a comparison function for f$(k), hence it follows from Lemma [2] that fs(k) > for 
k G (0, 1). □ 

Lemma 8. 1) Let u\ = ir, x G (0,1]; then sgn J\ = — sgn f zu (7T, k) = — sgn(2E(k) — K(k)), i. e., 
Ji < for k G (0, ko), Ji = if k = k , and J\ > for k G (fc , !)■ 
2) If u\ = ii, x = 0, i/ien Ji = 0. 

Proof. From a direct calculation it follows that 

J l (it 1 x,k) = -4x(l-k 2 x)f z (iT,k)h(k)i (23) 
/ z (7r,fc) = 2E(k)-K(k), 

f 3 (k) = E{k) + 2(k 2 - l)E(k)K(k) - {k 2 - l)K{k). 

Now the statement of item 1) of this lemma follows from Lemma [7] (fz(k) > for all k G (0, 1)) and 
the distribution of signs of the function 2E(k) — K(k) |llj (this function is positive for k < ko, equals 
zero if k = ko and is negative for k > ko). 

The statement of item 2) follows from formula (1231). □ 



Lemma 9. Let k G (0, 1), k / ko, u% = u\{k). Then a\[u\, k) < 0. 
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Proof. From a direct calculation it follows that if i±2 = u\(k), then 
ai = (eo + eiF(m) + e 2 F 2 (u 1 )) / (4 cos 3 ui) , 
eo = cos 2 u\\/l — k 2 sin 2 ui (l — A; 2 (l — cos 4 Ui)) , 

ei = — 2A; 2 cos ui sin tii (1 — (4 — 5 cos 2 ui + cos 4 u\) + (5 — cos 2 u\) k A sin 4 u\ — 2k 6 sin 6 u\), 

3 

e2 = sin 2 u\ (l — A; 2 sin 2 u%\ 5 (l — k 2 (l — cos . 

To estimate a sign of the function a\ notice first that u\ = u\ G (§, ^) , therefore cos 3 iti < 0. Further, 
we analyze a sign of the quadratic trinomial h(z) = eo + &\z + e2.z 2 . Its discriminant is equal to 

O 0/1 /I / O O \ 4 

D = ej — 4eoe2 = — 16A; sin tti cos u\ ( 1 — k sin txij , 
therefore D < for fc ^ fco- We have /i(z) > for A; 7^ ko because eo > 0, therefore a\ < 0. □ 

Lemma 10. 1) Lei A; G (0,1), Ui = u^(A;), x G [0,1); f/ien sgnJ x = sgn f zu (ir,k) = sgn(2E(k) — 
K(k)), i. e.j Ji > for k G (0, A; ), J\ = if k = k , and J\ < for k e (k ,l). 
2) If he (0,1), «i = u^jfe), x = 1, i/iera J x = 0. 

Proof. From equalities (fT7|) . ([TB]) , (fT6j) we get for m = uJ(As), i. e., if f zu {ui,k) = 0: 

d 4 = 0, 
d 2 = -do, 

Ji = do(l - x) = ai sin ul(k) (I - x). (24) 

It is shown in [11] that the function sinu3,(A;) = snp3,(A;) is negative for k G (0, A;o), is equal to 
zero if k = ko, and is positive for k G (Aio, 1). Thus for x G [0, 1) there holds the equality sgn J\ = 
— sgn ai-sgn(2E(k)—K(k)). To finish the proof of item 1) of the lemma we use Lemma[9j for u = u\{k), 
k ^ ko the function a\ is negative. Also, for k = ko, u = u],(ko = it) we have J\ = 2E{ko) — K(ko) = 0. 
Item 2) of the lemma follows from (JMJ- □ 

4.4. Global bounds of conjugate time in the subdomain C\. We prove estimate (jHJ) and get 

the upper bound for the first conjugate time in this subsection. 

Theorem 4. IfXed, then tl on - (A) > ^ A x( A )- 

Proof. Let A = (<p, A;, a = 1) G C\. 

1) Suppose k G (0, A;o)- It is required to show that ^ on j(A) > AK{k). 
1.1) Let sn((/?, k) 7^ 0. Consider the family of extremal trajectories 

q s (t) = Exp(A s ,i), te[0,tl\, 8€(0,ko), 

X s = (if s ,k s ,a = l) G Ci, 

A; 5 = s, ^ s = F(am(y?, fc),s), if = 42f(s). 

For any trajectory from this family 

X s = sin 2 ul = sn 2 t s = sn 2 ( V s + f ) = sn2 + 2i ^ S )> k °) = sn2 (V*. = 

= sin 2 (am (ip s , A; s )) = sin 2 (am(</?, A;)) = sn 2 (99, A;) / 0, 
'if 



am 



= a m(2K( S ),s)=7r, 



therefore from Lemma [8] J\(u\ , x s , A; s ) < 0. Namely the endpoint of a trajectory g s (if), s G (0, A;o), 
is not a conjugate point. According to Corollary CO there exists k G (0, A:), s. t. the trajectory q k (t), 
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t G (0, t\] does not contain conjugate points. We apply Theorem[3]to the family of the trajectories q s (t), 
s G [k,k] and see that the trajectory q k (t), t G (0, t\], has no conjugate points, i. e., t l co ^(X) > 4K(k). 
1.2) Let sn(c^, A;) = 0. Consider the family of trajectories 

q s = Exp(A s ,i), t e [0,t?], s€(0,e), 
A* = (c^ s ,£; s ,a = 1) G Ci, 



s, 93 



<^ + s, if=4X(s), 

where e > is a sufficiently small number, s. t. sn(c^ + s, k) 7^ for s G (0, e]. For the trajectories of 
this family we have 



x 



sn ip" + 



so 1 (up + s, k) ^ 0, sG(0,e], 



7T, 



therefore according to item 1.1) of this theorem, trajectories q s (t), t G (0, if], s G (0, e], have no 
conjugate points. Take any t\ G (0, 4K(k)). Since conjugate times are isolated from each other, 
it follows that there exists t% G (ti,4K(k)) that is not a conjugate time along the trajectory q°(t). 
Thus the instant ti is not a conjugate time for all trajectories of the family q s (t), s G [0,e]. Using 
Theorem^ we see that the trajectory q°(t), t G (0, £2] has no conjugate points. Therefore the instant 
t\ is not a conjugate time. Since t\ G (0,4K(k)), we obtain the required inequality iJ on j(A) > 4K(k). 
Note that the equality is attained in this case: from Lemma[8jit follows that J\(u\,x, k) = therefore 

4,nj(A)=4^(fc). 

2) Suppose k = ko- Take any t\ G (0,4if(fco)) and any t% G (ti,4K(ko)), which is not a conjugate 
time for the trajectory Exp(A,i). Applying Theorem to the family 

(f = Exp(A*,i), t€[0,tf], a€(-e,0), 

A s = (^ s ,A; s ,a = 1) G Ci, 

/c s = /c + s, ^ s = ¥>, t s i=h, 
where e > is sufficiently small as in item 1.2) of this proof, we see that ^onj( A ) ^ 4K(k ). According 
to LemmaEl ^ onj (A) > 4K(k ). 

3) Suppose k G (ko, 1). We claim that for any x G [0, 1] the set {(«i, k)\Ji(ui, x, k) = 0} is contained 
between the curves u\ = it and u\ = u\(k) in a neighborhood of (u\,k) = (71", &o); it can easily be 
checked that these curves are smooth and meet in the point (ui, k) = (it, ko) at the right angle. We 
have at this point the expansion: 

8E 3 (k 



Ji(ui,x, k) 



+x 



4E z (k )(iT-u l ) + x 
"k 



k (l 



■(k - k ) + 4E\k )(7T - Ul ) 



+ 



1 



k 2 



E\ko){k - ko) + 0((k - koY + (vr - Ul ) 2 ). 



Thus we get: 



x ^ 1 



1 



dJx 

~dk 



4E 2 (k )(x-l)^0, 
8E 3 (k ) 



kn 



ui=7v,k=ko 

Therefore for any x G [0, 1] the equation J\(u\, x, k) = defines a smooth curve in a neighborhood 
of the point (u±,k) = (w,ko). From Lemmas 151 and [TU1 it follows that for any x G [0,1], k G (0,1) 
the function J\(u\, x, k) equals zero on the interval u% G [it, u^(k)]. Therefore the curve {J\ = 0} is 
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contained between the curves {u± = 7r} and {u\ = u\(k)\ near the point (ui, k) = (ir, ko)- Hence for 
any x G [0, 1] there exists a neighborhood of the point (ui, k) = (n, ko), which satisfies the inequality 
J\{u\, x, k) 7^ for u\ < mm(-K,ul(k)) = ^max(^)- 

3.1) Let sn 2 ((/? + p\(k), k) ^ 1. For x = sn 2 (^p + ^MAS^j = gn 2 (c^ + p l z (k)), in a neighborhood O of 

the point {u\,k) = (n,ko) the function Ji(u\,x,k) does not vanish for k > ko, u\ < u\(k). Applying 
Theorem [3] to the family of trajectories 

q s {t) = Exp(A*,i), t€[0,tl], se[k ,k], 
X s = (p s ,k s ,a = 1) G Ci, 
k s = s, ip s = F(&m(ip +p l z (k),k),s), 
4 = 2F(u 1 z (k),s), 

we see that the trajectory q(t) = q k (t) = Exp(A,i), A = X k = (ip, k, a = 1), t G [0,t\, t\ = t k , has no 
conjugate points. 

Finally, applying Theorem [3] to the family 

q s (t) = Ex P (A s ,t), te[0,tf], se[fc,4 

A* = (v s ,k s ,a = 1) G Ci, 

k s = s, tp s = F(axa(tp+pl(k),k),s), 

t{ = 2pl(s), 

we see that the trajectory q k (t) = Exp(A,i), t G [0, ijyiAxM] 113,8 no conjugate points, Q. E. D. 

3.2) In the case sn 2 (c/? + p\{k), k) = 1, the proof of t\ on ^{\) > 2p\{k) is obtained as in item 1.2). 
The theorem is completely proved. □ 

Remark 1 . The lower bound in estimate of conjugate time (|8|) is attained. From Lemmas El [10] and 
Theorem 2] we get: 

k e (0, ko), simp = 0, a = 1 t^-(A) = 4JT(A;) = 4axW. 

fc = fc , a = i => 4 n .(A) = 4K"(fco) = 4ax( a ) ; 

fc € {ko, 1), sin 2 (^ + p*(A:)) = 1, a = 1 =► 4 nj (A) = 2p\{k) = iJjAxW- 

In addition to the lower bound from Theorem H] we get the upper bound of the first conjugate time 
in terms of the second Maxwell time 

4ax(A) = max(2pi, 4K)/a, A G C v 

For a = 1 we obtain: 

fcG(0,& ) => t 2 MAX (X) = 2p 1 z (k), 
k = k => tMAx(A) = 2pi(*o) = 4*r(fco), 
fc€(fco,l) => 4ax(A) = 4iT(fc). 
Proposition 2. If X £ C\, then t l con ^X) < *m AX (A). 

Proof. Let A = (93, fc, a) G C\. Also, in the proof of the lower bound of conjugate time we can assume 
that a = 1. 

If k G (0, ko), then for any x G [0, 1] we have Ji(n,x,k) < (see Lemma[8|) and Ji(uj.(/c), x, fc) > 
(see Lemma flOj) . i. e., the function J\{u\,x,k) changes sign in the segment u\ G [tt, u\{k)]. Therefore 
the corresponding segment t G [4K(k), 2p\{k)\ = [*max(^)' ^max(^)] contains the first conjugate time. 

If k = ko, then for any x G [0, 1] we get Ji[tt,x, k ) = (see Lemma|S| thus tl on ^(X) = AK(k ) = 

^MAx(^) = ^MAx(^)- 
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Finally, if k > ko, then for any x G [0,1] J{u\{k),x,k) < (see Lemma [TO)) . Ji(-K,x,k) > (see 
Lemma©, therefore ^ onj (A) G [2p\{k), 4K (fc)] = [i^xW, ^ AX (A)]. □ 

Remark 2. One should not think that the segment [^max(^)' ^max(^)] contains exactly one conjugate 
time. Computational experiments in the system Mathematica show that for if = and k G (0, 999; 1) 
this segment contains two conjugate times. 



5. Estimate of conjugate time for A G C 2 

The aim of this section is to prove estimate ([8]) in the domain C2 for a = 1: iconj(^) — 2Kk, A G 
Using parameterization of extremal trajectories [5] for A = (ip, k, a) G C2, as well as in the domain C\ 
we get the expression of the Jacobian J = g^'^'a) for a = 1: 

J = R ■ i/i, 

R= 32 _^ 

k(l — k 2 )(l — k 2 sin 2 u\ sin 2 U2) 2 ' 
J\ = do + d2 sin 2 U2 + sin U2, 
Ml = am(p, fe), U2 = am(r, k), 
t _2ip + t 

P ~Zk' T ~^k~' 



d = -sinuicostii(4E ,3 (ni)sin(2ni) -4(1 - k 2 ) cos(2ui)\/l - k 2 sin 2 mF 2 (m) + 

+(8 - 8/c 2 + A; 4 + fc 2 (2 - k 2 ) cob(2«i)) sin(2ui)F(«i) + 2(2 - 3A; 2 + A: 4 ) sin(2ui)F 3 (m) + 
+2/c 2 a/1 - A; 2 sin 2 u[ sin 2 (2m) + 2F 2 (m)(6cos(2m) \A - k 2 sin 2 u[ - (2 - A; 2 )sin(2ui) x 
xF(«i)) - E(«i)(4(2 - A: 2 ) cos(2ui)\A - & 2 sin 2 mF(m) + 2(4 - 2A; 2 + 
+3k 2 cob(2«i)) sin 2m + 4(1 - k 2 ) sin(2Mi)F 2 (m), 



c?2 = —2k 2 (l — k 2 ) cos U\ sin 3 ni\/l — A; 2 sin 2 uiF 2 (m) — 2/c 4 cos 3 m sin 3 u\ yl — k 2 sin 2 u\ — 
-2(1 - A: 2 sin 4 Ul )E 3 ( Ul ) - (2 - 3/c 2 + A; 4 )(l - A; 2 sin 4 Ul )F 3 ( Ul ) + 
+E 2 ( Ul ){6k 2 cos u\ sin 3 m \/l - A; 2 sin 2 m + (2 - A: 2 )(l - A: 2 sin 4 m)F(m)) + 
+F(ui)(A; 2 cos 2 ui sin 2 m(4 - 3A; 2 + 3A; 2 cos (2m)) - 

-2A: 2 (2 - A: 2 ) cosm sin 3 my/1 - k 2 sin 2 mF(m) + 2(1 - A; 2 )(l - A: 2 sin 4 u 1 )F 2 {u 1 )) - 

--A; 2 (8 - 8A; 2 + A: 4 + k 2 {2 - k 2 ) cos(2m)) sin 2 (2m)F(m), 
8 

d 4 = 2(1 - A: 2 sin 2 m)F 3 (m) - F 2 (m)(3A; 2 cos m sin m Vl - A: 2 sin 2 m + (2 - A; 2 )(l - 
-A: 2 sin 2 m)F(m) + t(1 - A; 2 )F 2 (m)(2(2 - A; 2 )(2 - A: 2 + A: 2 cos(2«i))F(«i) + 

+2A; 2 \/l — A; 2 sin 2 m sin 2m) + — F'('Ui)(4fc 4 cos 2 m sin 2 u\ — 

-8(1 - k 2 )(l - k 2 sin 2 m)F 2 (m) + 2A: 2 (2 - A: 2 ) sin(2m)\/l - ^ 2 sin 2 mF(m))- 



14 



5.1. Conjugate time as k — > 0. Asymptotics of the function J\ as k — > has the form: 

k 8 

Ji(ui,x,k) = ——^jQ(ui,x) + o(k 8 ), x = sm 2 u 2 , (25) 
J?(u!,x) = dg(m) + c^(m)a + d^(ui)x 2 , 

dg(wi) = - cos u\ sinni ((— 48uf — 3) cos(2«i) + 3 cos(6ui) + (42«i — 64uf ) sin(2 / ui) + 2u\ sin(6ui)) , 
8 

^2(^1) = -^4(^1) = -(sin(4ui) - Au\)(Au\ + sin(4ui)iii + cos(4ni) - 1). 
First we prove several auxiliary lemmas which give an estimate of the functions d%. 
Lemma 11. The function fi(u\) = 8u\ + Au\ cos(4ui) — 3sinui > on the interval u\ G (0, 5). 



Proof. We show that the function g{u\) = 2 + cos(4ui) is a comparison function for f\(u\) on the 

cpansion fi(u\) = ^§-u 5 + < 
/i(m) V = 16sin 4 (2m) 
g{ui) ) (2 + cos(4ui)) 2 



interval u\ G (0, ^). The inequality f\(u\) ^ follows from the expansion fi(u\) = ^§-u 5 + o(u 5 ). 

xt i / \ rv r „. „ , ,. . f fi(ui)\ 16sin 4 (2«i) 

Note that > tor any u\. finally we get the equalities I — - — r- = — — — r^r and 



g(ui) 45 

So g(ui) is a comparison function for fi(ui) thus it follows from Lemma [2] that > for 

ui €((),§). □ 

Lemma 12. T/ie function f2(u±) = —1 + 4u 2 + cos(4rti) +ui sin(4ui) > on i/ie interval u\ G (0, ^). 

Proof. Let us show that the function <?(ui) = 4ui + sin(4ui) is a comparison function for f2{ui) on 
the interval ui G (0, ^). The inequality ^(^l) ^ follows from the expansion ^(^l) = ipf^ 6 + o(u 6 ). 
If tii > 0, then «i + sinni > 0, therefore (7(1*1) > for u\ G (0, -|). Finally we have the equalities 
h{ux)\ _ (-4m+sin(4m)) 2 mJ / 2 (m) _ 16^ , 5 



, , ,, and — ; — r- = — iii + o(ui). So g(iii) is a comparison function for 

S(ui) y (4«i +sm(4 Ul )) 2 g( Ul ) 45 

f2{u\) thus it follows from Lemma [2] that f 2(1*1) > for ui G (0, f ). □ 

Lemma 13. 7/ui G (0,^), i/ien the function d%{u\) = | cosui sinui((— 48uf — 3) cos(2ui) + (42«i — 
- QAuf) sin(2«i) + 3 cos(6ui) + 2u\ sin(6ui)) > 0. 

Proof. We now prove that the function g(u\) = 4«i+sin(4ui) is a comparison function for oIq(ui) on the 
interval u\ G (0, |). The inequality d^(u±) ^ follows from the expansion dg(wi) = |^|| u 11 + o(u n ). 
If u\ G (0, ~), then > 0. In the equation 

' /l(ui)/2(«l), 



<?(ui) / 8uf cos 2 tti sin 2 ni 
we note that f\(u\) = 8«i + 4ui cos(4u!) — 3sin(4ui) > on the interval Ui G (0, |) (see Lemma [TT 



and f2(ui) = — 1 + 4-u 2 + cos(4ui) + uisin(4ui) > on the interval u\ G (0, |) (see Lemma fT2j) , 
Meanwhile -j^j = 37§f' u i + °( u i)- ^° 9( u i) ^ s a comparison function for ^(m), therefore it follows 
from Lemma [2] that d[](ni) > for m G (0, f ). □ 

Now we estimate the function J®. 

Lemma 14. For any u\ G (0, x G [0, 1] the inequality J®(ui,x) > holds. 

Proof. It follows from Lemma [H that J? («i,0) = J? 1) = dg(«i) > for all m G (0, f). Further, 
it follows from Lemma [TZI that d\(u\) = |(sin(4ni) — 4«i)/2(«i) < for U2 G (0, ^). Therefore the 
statement of this lemma follows from Lemma [3l □ 
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Proposition 3. There exists k G (0, 1), s. t. for any k G (0, k), U\ G (0, x G [0, 1], i/ie inequality 
Ji(ui,x,k) > holds. 

Proof. This proposition is proved in exactly the same way as Proposition [U with the use of Lemma [LA 
expansions (125p and the following expansions: 



J x = -J-fcX 1 + +o(k 8 u\ 1 ) + ^=k % u\x + o(A: 8 n?x) - T^A; 8 n?x 2 + o(k*u\x 2 ), k 2 + u 2 ^ 0, 
Ji = |rfe^ 8 (f -u 1 )+o(k s (^ - Ul )) + ^k 8 x +o(k 8 x)- 

-^k s x 2 + (k 8 x 2 ), k 2 +(^- ul y^o. 



□ 



From Proposition [3] we get the following statement in the variables (t,(p,k). 



Corollary 2. There exists k G (0,1), s. t. for any k G (0, k), ip G R, the trajectory Exp(A,t), 
A = ((p, k,a) G C2, t G (0,i^ AX (A)), does not contain conjugate points. 

5.2. Conjugate time for t = ij^AX - ^he instant of time ^max(^) = corresponds to the value 
of the variable ui = | . We have 

TT \ TT 7T 



Ji (^-, x, fcj = d%x + dl x 2 , (26) 

4 = -dl = ^/l-k 2 g z (K, k) U(k), (27) 
g z (p,k) = ((k 2 - 2)p + 2E(p))dnp- k 2 snpcnp, (28) 
U{k) = E 2 {k)-{l-k 2 )K{k). (29) 

In the paper [TT| it was proved that g z (p, k) < for any p > 0, A; G (0, 1); therefore g z (K, k) < 0. 

Lemma 15. TTie function fi{k) = E 2 (k) + (A; 2 — l)/r(/c) > on i/ie interval k G (0, 1). 

Proof. We show that the function g(/c) = 1 — k 2 is a comparison function for f^(k) on the interval 
k G (0,1). The inequality fn{k) ^ follows from the expansion U(k) = §2 A; 4 + o(k 4 ). Note that 
g(k) > for /c G (0, 1). Finally we have the equalities 

[U(k)\ = 2(^(fc) + (fc 2 - l)K(k)) 2 
V 9(k) J k(k* - l)^ 

fi{k) TT 2 

and — j-^y = + o(k 4 ). So <?(&) is a comparison function for f^{k) thus it follows from Lemma [2] 

that / 4 (fc) > for k G (0,1). □ 

Lemma 16. 1) If k £ (0, 1), uj = |, 1 € (0, 1), ^en Ji > 0. 
2) If k £ (0,1), ui = § , x G {0,1}, then J x = 0. 

Proof. It follows from formula ([26]) - ([29]) . inequality g z (K,k) < 0, and Lemma [15j □ 

5.3. Global bounds of conjugate time. 
Theorem 5. //AG C 2 , then 4 nj (A) > ^ A x( A )- 
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Proof. This theorem is proved in exactly the same way as Theorem U] based on homotopy invariance 
of index of the second variation (the number of conjugate points), see Theorem [3j The last theorem 
is applied to the family of extremal trajectries 

q s (t) = Exp(A s ,t), i€[0,if], s€[k,k], 
X s = (ip s ,k s ,a = l) G C a , 

<p* = sF (am ^ + ^ /2 , fcV - if/2, 

k s = s, 4 = 2K(s)s, k£{0,k). 

□ 

Remark 3. It follows from Lemma[l6\that the lower bound from Theorem\5\is attained: if ip = Kkn, 
neZ, then t\ on] {\) = 2Kk, A = (<p, k, a = 1) G C 2 . 

Remark 4. Using the homotopy invariance of the index of the second variation, we can prove the 
upper bound of conjugate time: 

iconj(^) < *MAxW: ^ e ^2, 

4ax(A) = 4fcK, A G C 2 . 

Ao£e t/iat the segment [£max(^)> *max(A)] contains exactly two conjugate times (with account of mul- 
tiplicity). 

6. Estimate of conjugate time for A G C3 

Theorem 6. If X £ C3, then the extremal trajectory Exp(A, t), t G (0,+oo) ; does not contain conju- 
gate points. 

Proof. Let A = (p, k = l,a = 1) G C3 and ti > 0. We show that the trajectory Exp(A, t) , t G (0, t{\, 
has no conjugate points. Choose a time i 2 > ti that is not a conjugate time. There exists k\ G (0, 1), 
s. t. k\K(ki) > 2t2- According to Theorem[5j all trajectories 

q s (t)=Exp(X s ,t), i€(0,tj], s€[fci,l), 
A s = (/^> = i)ec 2 , 

^ = ^ /c s = s, i? = -K(s)s, 

have no conjugate points. Applying Theorem[3]to the family of trajectories q s (t), t G (0, £2], s G 1], 
we conclude that the trajectory = Exp(A, i), i G (0,i 2 ], does not contain conjugate points. □ 

7. Estimate of conjugate time for A G uJ =4 Cj 

If A G uJ =4 Cj, then conjugate time (and cut time) can be located by projecting the original prob- 
lem {I])-© into simpler problems of a lower dimension using the following proposition. 

Proposition 4. Let us consider two optimal control problems: 

q l = f(q\u), q i G M i , u G U, 

q i (0)=ql q i (t 1 )=q\, 
h 

J = J (p(u) dt — V min, 


i = 1,2. 
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Suppose that there exists a smooth map G : M 1 — > M 2 , s. t. ifq l {t) is the trajectory of the first system 
corresponding to a control u(t), then q 2 (t) = G(q 1 (t)) is the trajectory of the second system with the 
same control. 

Further assume that and q 2 (t) are such trajectories. If q 2 (t) is locally (globally) optimal for 

the second problem, then is locally (globally) optimal for the first problem. 

Proof. Assume the converse. Suppose q 2 (t) is optimal and is not optimal. Then for the first 

problem there exists a trajectory q (t), s. t. value of the functional J for this trajectory is less than 
for So value of J is less on the trajectory q 2 (t) = than on q 2 (t). This contradiction 

proves the proposition. □ 



Riemannian problem in the Euclidean plane 10 2 



In the case A G C4 U C5 U C7, the sub- Riemannian problem on the Engel group is projected on the 

2 . 

x,y 

G : -> {x, V, z, v) i-> 0, y), 

x = ui, y = u 2 , (x,y)(0) = (x ,yo), (x,y)(h) = 



ti 



If A G C4 U C5 U Cj, Exp(A, t) = (xt,yt, zt,v t ), then (xt,yt) is a straight line that is globally optimal 

l 2 x , y for t G [0,+oo). Therefore i cut (A) = 4 nj ( 



in the Riemannian problem on for i G [0, +00). Therefore i C ut(A) = tl-(X) = +00 = %Ax(^) 



for A G C 4 U C 5 U Oj. 

If A G Cq, then the sub- Riemannian problem on the Engel group is projected on the sub-Riemannian 
problem on the Heisenberg group M%. y : 

G ■ K,y,z,v -> R l,y,zi 0> V, z , v ) ^ 0> V, z ), 

y x 

x = ui, y = u 2 , z = --ui + -u 2 , 

(x, y, z)(0) = (x , y , z ), (x, y, z)(t x ) = Z\), 

I = / a / Vn ~\~ Mn (it — y min . 



For A = (9, c, a = 0) G Cq, Exp(A, t) = (x t , yt, zt,vt), the curve (xt, yt, Zt) is globally and locally optimal 

, . / cos (ci + 6) — cos sin(ci) — sin( 

tor t G 10, t4|, 1. e., up to the hrst turn of the circle (xt,yt) = , 

|c| \ c c 

It follows from Proposition H] that tl on ^(X) > t cut (X) > ^ = ^max(^) f° r ^ ^ C§. By Theorem [TJ we 



6- 



have ^ onj (A) > i cut (A) < g = ^ AX (A) for A G Q 

Remark 5. Passing to the limit a — > 0, k — > 0, it can be shown that for A = (9, c, a) G Cq, 6 = a = 0, 
equality i,l on j(A) = ]tf = ^max(^) holds. But for A G Cq this equality does not hold in the general case. 

Finally we summarize the results of this section in the following statement. 

Theorem 7. If X G C4 U C5 U C 7 , then ^(A) = *cut(A) = +00 = t^^X). If A G C 6 , then 

*conj(^) - *cut(A) = i M AxM- 

8. Conclusion 

Theorem [2] follows from Theorems IH El El El 

Using the estimate of cut time obtained by the work [5] (Theorem [TJ and the estimate of conjugate 
time proved in this work (Theorem^, we can get the description of global structure of the exponential 
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map in sub-Riemannian problem on the Engel group. So we can reduce this problem to solving the 
system of algebraic equations. This will be the subject of another paper. 

The method for estimating a conjugate time used in this paper was successfully applied earlier to 
Euler's elastic problem [7] and sub-Riemannian problem on the group of rototranslations [12]. There 
is no doubt that this method is also valid for nilpotent sub-Riemannian problem with the growth 
vector (2,3,5) [EHD]. 

The method can be used for other invariant sub-Riemannian problems on Lie groups of low- 
dimensional integrable in non-elementary functions. 

The first natural step in this direction is investigation of invariant sub-Riemannian problem on 3D 
Lie groups which are classified by A.A.Agrachev and D.Barilari [2]. 
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